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ABSTRACT 

Stimulated by recent observations of a giant radio-to-X-ray synchrotron flare from 
HST-1, the innermost knot of the M 87 jet, as well as by a detection of a very high 
energy 7-ray emission from M 87, we investigated the dynamics and multiwave- 
length emission of the HST-1 region. We study thermal pressure of the hot interstel- 
lar medium in M 87 and argue for a presence of a gaseous condensation in its cen- 
tral parts. We postulate that this additional feature is linked to the observed central 
stellar cusp of the elliptical host. Interaction of the jet with such a feature is likely 
to result in formation of a stationary converging/diverging reconfinement/reflected 
shock structure in the innermost parts of the M 87 jet. We show that for a realistic set 
of the outflow parameters, a stationary and a flaring part of the HST-1 knot located 
^ 100 pc away from the active center can be associated with the decelerated portion 
of the jet matter placed immediately downstream of the point where the reconfine- 
ment shock reaches the jet axis. We discuss a possible scenario explaining a broad- 
band brightening of the HST-1 region related to the variable activity of the central 
core. In particular, we show that assuming a previous epoch of the high central black 
hole activity resulting in ejection of excess particles and photons down along the jet, 
one may first expect a high-energy flare of HST-1 due to inverse-Comptonisation of 
the nuclear radiation, followed after a few years by an increase in the synchrotron 
continuum of this region. The synchrotron flare itself could be accompanied by a 
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subsequent inverse-Compton brightening due to upscattering of the ambient (mostly 
starUght) photons. If this is the case, then the recently observed order-of-magnitude 
increase in the knot luminosity in all spectral bands could be regarded as an unusual 
echo of the order-of-magnitude outburst that had happened previously (and could be 
eventually observed some ^ 40 years ago) in the highly relativistic active core of 
the M 87 radio galaxy. We show that very high energy 7-ray fluxes expected in a 
framework of the proposed scenario are consistent with the observed ones. 

Key words: radiation mechanisms: non-thermal — shock waves — galaxies: active 
— galaxies: jets — galaxies: individual (M 87) 



1 INTRODUCTION 

A kiloparsec-scale jet in M 87, t he giant ellip tical radio galaxy in Virgo cluster — the very first 
extragalactic jet ever discovered JcurtiJIlQlsl — provides us with an exceptional laboratory for 
studying physics of relativistic coUimated outflows. This is because its proximity (16 Mpc, 1 arc- 
sec = 78 pc) allows for observations of the jet and of ambient medium at different frequencies, 
with a very high spatial resolution. This jet has in fact been studied in great detail in all wave- 



lengths. One of the most remark able features of this jet is the inner HST 



about 60 pc from the active core ( Biretta et al.ll999 ; Perlman et ahLoOS ; Harris et ahLoOS 



knot region, observed at 



Complexity of this innermost part of the outflow, which consists of stationary and superluminal 
sub-components, as well as its uniquely variable broad-band emission, calls for an explanation. 
Here we investigate the properties of HST-1 knot and present an attempt to provide such an expla- 
nation. 

Below we summarize observational constraints on the physics of the M 87 jet in general. 
Next, in section 2, we investigate properties of the ambient medium necessary for understanding 
the dynamics of the jet. We argue that the enhanced thermal pressure connected with a stellar 
cusp in the innermost parts of the host galaxy is likely to form a stationary converging/diverging 
reconfinement/reflected shock structure in the jet. In section 3, we conclude that indeed the sta- 
tionary and flaring upstream edge of HST-1 knot can be associated with the decelerated portion 
of the jet matter placed immediately downstream from the point where the reconfinement shock 
reaches the jet axis. The presented interpretation gives us a framework in which one can model 
the broad-band emission of this part of the jet, and its high energy 7-ray emission in particular. In 
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section 4, we present an evaluation of the radiation fields at the position of HST-1 knot. We find 
that the energy density of the starlight and of the nuclear emission can be comparable to the en- 
ergy density of the equipartition magnetic field in this jet region. This means that the high energy 
7-ray inverse-Compton emission of HST-1 knot's flaring point should be expected at the (very 
roughly) similar level as its observed radio-to-X-ray sync hrotron emission. Since the latter one 



20061) ■ the expected TeV 



has recently increased significantly up to ~ 10"^^ erg s ^ dHarris et al. 
flux fro m HST- 1 knot should then be promisingly c omparable to the one detected from the M 87 



system dAharonian etal .112003 ; 



Beilicke et alJl2005|) . This issue is investigated further in section 5. 



There we present a possible scenario relating variable emission of the HST-1 knot/reconfinement 
nozzle with a modulated activity of the relativistic central core. In particular, we show that assum- 
ing the previous epoch of the high nuclear activity of a central black hole, resulting in ejection 
of excess particles and photons down along the jet, one may expect first a high-energy flare of 
HST- 1 due to inverse-Compton scattering of the nuclear radiation, followed a few years later by 
an increase in the synchrotron continuum emission of this region. Interestingly, the predicted 7-ray 
fluxes (assuming energy equipartition between radiating electrons and the jet magnetic field) are 
consistent with the observed ones. Final summary and conclusions from our study are presented 
in the last section 6. 



1.1 M 87 Jet 



Junor et al. 



(119991) show a presence 



Very Long Baseline Interferometry observations reported by 
of a very broad radio-emitting limb-brightened outflow close (~ 10~^ parsec) to the M 87 center, 
characterized by an opening angle of $obs ~ 60°. This outflow experiences a strong coUimation at 
the projected distance equivalent to ~ 100 (Schwarzschild radii) from the central supermassive 
black hole (hereafter 'SMBH'), i.e., with the appropriate conversion = 3.85 yuarcsec = 0.003 pc, 
at the distance of ~ 0.4 mas ^ 0.03 pc. The coUimation continues out to ^ 10 pc from the center, 
where the jet adopts an opening angle $obs ^10° that remains roughly stable further away from 
the core. A detection of synchrotron self-absorption features in the radio spectrum of the M 87 
nuclear regio n allowed the placem ent of an upper limit on the jet magnetic field B < 0.2 G at 



0.06 pc (R eynolds et al 



the radio spectra. 



1996r). Re centlv. by modeling a turn-over frequency along the jet in 



Dodson et al 



(I2005h found B ~ 0.01 - 0.1 G for r < 3 mas ^ 0.25 pc, and 
B < 0.01 G further away along the jet. 

The bright components of the radio jet placed at r ~ 0.1 — 5 pc from the jet base are char- 
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acterized by no, or some sub-relativistic proper motions, /3app ^ 0.04 (IJunor & Biretta 



1995 



Dodson et al.ll2005f) . Further out along the jet, the knots detected at several pc from the core (e.g., 



the knot L placed at r ~ 0.16" ~ 1 2.5 pc) are also s ub-luminal. 



apparent velocities f3., 



app 



0.3 - 0.6 dReid et al 



1989 : 



loweyer with slightly larger 



Biretta et al. 



19991) . Surprisingly, several 



knots placed even further out (0.8" — 6.3") occurred to be highly su perluminal, with 0^^ ^ reach 



Biretta et al. 



( 19991) showed 



ing 6. In particular, Hubble Space Telescope observations reported by 
that the unresolved stationary feature upstream of the HST-1 knot (at 0.8" ~ 62 pc) seems to 
emit various components down the jet, both slow and fast, with the maximum apparent velocity of 



app 



5 — 6. Also, all the components of knot D (2.7" — 4" ^ 210 — 312 pc) are superluminal. 



with /3app ~ 2.5 — 5. Finally, knot E placed at r ~ 6" ^ 500 pc from the core is characterized 
by a relatively high velocity of /5app ~ 4. All the measured superluminal features, if interpreted as 
moving blobs, suggest the bulk Lorentz facto r for the 1" — 6" po rtion of the jet larger than F ^ 6, 
and the jet viewing angle less than 6 ^ 20° iBiretta et al.lll999|) . (With 9 ^ 20° all the projected 
distances along the jet cited in this section should be multiplied by a factor of 3.) 

A stationary feature placed at the upstream edge of HST-1 knot has been flaring in the op- 
tical and X-rays since 2002. The results of Very Large Array ^ Cha ndra X-ray Observato ry and 



Hubble monitoring programs presented by 



Harris et al. 



O003) and 



Perlman et al. 



( 20031) estab- 



lished month-to-year variability of its radio-to-X-ray synchrotron continuum, with a comparable 
amplitude over the entire broad waveband. The HST-1 knot is unresolved by Hubble, indicating its 
spatial dimension smaller than R ^ 0.02" ^1.5 pc. The equipartition magnetic field at the posi- 
tion of this knot, when evaluated at the quiescence state and neglecting corrections due to expected 
relativistic bulk velocity of the radiating plasma, is of the order of i?cq ~ 10~^ G. Magnetic field 
lines thereby are pred ominantly perpendicular to the jet axis, as suggested by polarization studies 



(IPerlman et al 



2003h . The degree of the linear polarization decreases from 0.68 upstream of the 



HST-1 flaring region (0.72" ^ 56 pc), to 0.46 at the position of the flux maximum (0.8" ~ 62 pc), 
and then to 0.23 downstream of it (0.92" ^ 72 pc). The most recent data sho w that until the yea r 



2005 the X-ray emission of HST-1 knot increased by as much as a factor of 50 (IHarris et al 



2006h 



At the distance r ~ 12.2" 950 pc away from the center the jet brightens significantly, form- 
ing a prominent knot (knot A) followed by the subsequent knots (B and C), to disa ppear at ap 



proxi mately r ~ 2 kpc into an amorphous radio lobe visible at low radio frequencies (lOwen et al 



20001) . VLA studi es indicate sublum inal apparent velocities of these outer jet components, with 



/Sapp < 0.5 — 0.6 (IBiretta et al.lll995l). The kpc-scale let, when observed in radio, exhibits filamen- 



tary limb-brightened morphology (lOwen et al 



19891) . Both its optical and X-ray emissions are syn- 
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chrotron in origin, similar to the inner parts of the ou t flow (iBiretta et al.lll991 



Sparks et al 



Perlman et al 



2m 



Marshall etal. 



Meisenheimer et al 



Wilson & Yang 



mm 



All along the jet the radio-to-optical power-law slope is 0.65 < q;r-o ^ 0.8, while the optical- 
to-X-ray one 1.0 < ao_y < 1.9 with the exce ption of HST-1 knot, for which ao-x ~ 0.8 — 1.0 



( Perlman & Wilsonl2005 



I 



Waters & ZepJ2005l) . This indicates a general 'broken po\ yer-law' char- 



Stawarz et al 



acter o f the broad-band synchrotron spectrum along the M 87 jet. As discussed in 
(120051) . the kpc-scale jet's magnetic field is not likely to be lower than B^q ~ 300 /xG. 

HEGRA Cherenkov Telescopes System detected the M 87 em ission with the photon flux of 
F^{> 0.73 TeV) 0.96 x 10"^^ cm"^ s"^ ( Aharonian et akkoOB ). Assuming a spectral index for 
the observed emission = 2, this corresponds to the isotropic of luminosity (0.73 TeV) ~ 
6.9 X lO''" erg s^^. The obseryations were taken in the period 1998-99, when the HST-1 flaring 
region was in its quiescence epoch. Different scenarios were proposed to account for the detected 
TeV signal, inc 



uding various yersions of moc 


eling M 87 active nuc 


eus ('i 


nisaligned' and 'struc- 


Bai & Lee 


2001; 


Reimer et al. 


2004; 


Ghisellini et al. 


2005; 


Georganopoulos et al. 



or of the kpc-scale jet (of its brightest knot A in particular; 



tured' BL Lac; 

2OO5L but also a high energy emission of the M 87 host galaxy (iPfrornmer & Ensslin 



200 



Stawarz etal. 



J 



120021). The eyidence 



for a year-timescale variability estab 



(LeBohec etal 



Beilicke et al. 



ished by the subsequent Whipple and H.E.S.S. observations 



2005|) . indicating a likely decrease of the TeV signal from 
M 87 from 1999 tiU 2004 by about an order of magnitude (down to L^(0.73 TeV) ^ 0.54 x 10^° 
erg s^^), excludes the later two po ssibilities , imposing however interesting constraints on the kpc- 



scale jet magnetic field intensity (IStawarz et al 



20051) . At the same time, synchrotron radio-to- 



X-ray emission o f the HST-1 flaring region increased by more than an order of magnitude (see 
Harris etaklllOOfil). 



2 HOST GALAXY EMISSION PROFILES 



Chandra studies presented by lYoung et al.l (120021) demonstrate that the X-ray surface brightness 
of Virgo A cluster (centered at the position of M 87 radio galaxy) is of the modified King profile 
Ex(r) oc [1 + (r/rx)^]"'^^"*""'^, with the parameter P = 0.4 and the critical radius tk ~ 18". 



This implies a density profile of the X-ray e mitting hot gas pcir) oc [1 + (r/rx) ] 



21-3/3/2 



I.e. 



oc r 



^1.2 



for r > rx (see, e.g.. 



Sarazin 



19861) . Both the temperature and the abundance of this 



gas decrease smoothly toward the cluster center, rea ching at r < 60" values of IcTq < 1.5 keV 



and Z < 0.5 Zq, respectively (IBohringer et al. 



20011) . It is not clear if the abundance decrease is 
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real, or only apparent, caused by resonant line scattering (but see 



Gastaldello & Molendi 



2002) 



However, even with this ambiguity, one can conclude from the X-ray observations that the number 
density of hot thermal electrons in the M 87 host gala xy decreases from ~ 0.15 cm~^ at r ^ 30 ^^ 



to ~ 0.03 cm ^ at r ~ 100" from the active core (lYoung et al. 



12002 : 



Pi Matteo et al. 



2003). 



All of these constraints are in agreement with a general finding that the central electron number 
density in g iant ellipticals is typically ~ 0.1 cm^^, and declines as oc r^^^^ with the distance from 
their cores ( Mathews & BrighentJ2003 ). Unfortunately, even with the excellent spatial resolution 
of Chandra the thermal gas X-ray emission profile cannot be probed precisely in the innermost 



portio ns of M 87, r < 10", because of numerous X-ray emission features there (see 



Feng et al 



20041) . As argued below, one can instead use the optical observations to constrain the parameters 



of this gaseous medium. 



Optical observations of M 87 reported bv lYoung et al. 



( 1978b indicate that the modified isother- 



mal sphere model (usually applied to elliptical galaxies) is inconsistent with the observed starlight 
emission profile for the projected ra dii of r < 10". In p articular, they showed a presence of a cen- 
tral luminosity excess, explained by lYoung et al.l (|l978|) in terms of a dynamical effect of a SMBH 



place d at the center of the galaxy on its stellar neighborhood. Hubble observations (ILauer et al 



1992h confirmed the presence of this additional stellar component in agreement with the interpreta- 
tion involving M ^ 3 x 10^ Mq SMBH strongly bounding nearby stars and creating a central stel- 



lar cusp with an increased ste llar velocity dispersion (in this context see also 



mi 



Dressier & Richstone 



Macchetto et al.ll997|) . The observed optical (starlight) surface brightness profile of the host 



galaxy is therefore described by So(r) oc with b = 0.25 for r < 3" and 6 = 1.3 for r > 10", 
and composed of two separate components, namely a central cusp with the luminosity density 
> 10^ Lq pc~^ (in the / filter for r < 0.1"), and a modified King profile with the curvature radius 
rc ~ 7" and a tidal radius tt > 10^ re- The later component is con sistent with a general property 
of the elliptical gala xies, na r nely lo g; rx/rc ~ 2.2 (jSilva et al 

]H). 



As discussed in 



Young (|l980l) . in the case of the adiabatic growth of a central SMBH (i.e.. 



the growth at a rate slower than the dynamical time scale of the stellar cluster but faster than the 
e), the expected den sity profile of stars follows a power law, psir) oc r^"- with 



relaxation time sea 
a = 1.5 (see also 



van der Marel 



19991) . Hence, the expected starlight brightness profile for the 



adiabatic stellar cusp is Snfr ) oc r x 



M 87. However, as noted by 



V) oc r i.e. much steeper than the one obse r ved in 



Uaueretal 



Young et al 



( 1992|) . the behavior discussed by 
is in fact an asymptotic one, holding for r — > 0, while in the outer regions of a cusp an expected 
profile should be flatter, similar to oc r~°^^ observed in M 87. This central stellar cusp profile joins 
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smoothly with the galactic starlight profile So('^) oc r at further distances from t he core, im 



plying a stellar density dependence of psir) oc r Note that Hubble observations (|Lauer et al 



plyinj 


J a stellar density 


199.4 


Ipaheretal. 


1997i 



of elliptical galaxies show that they never possess a homogeneous surface 
brightness distribution So(r) oc const expected in the case of a pure King-like profile, but they 
can be divided into two classes: (i) the 'core type' galaxies with the brightness profile described 
by a broken power law with So(r) oc r^^ and b ^ 0.3 within some critical (break) radius r < 
(in M 87 case re ~ 3"), and (ii) the 'power-law type' galaxies characterized by a sin gle value of 



deRuiter et al 



( 2n05n . radio loud 



b ^ 0.5 within the whole central region. As noted recently by 
galaxies are always of the 'core type', although not every 'core type' galaxy is radio loud. Also, 
the break radius is proportional to the galactic luminosity, roughly re/kpc ~ Ly/lO^^ergs"^. 
One can therefore suspect the presence of an additional thermal pressure component within the 
host galaxies of radio-loud AGNs, 'matching' the central stellar cusps. 



Trinchieri et al 



mm. 



Our main assumption follows from the observational fact reported by 
that the optical and X-ray surface brightness profiles for three bright Virgo A ellipticals are al- 
most identical, i.e. that So(r) oc Sx('")- Since the starlight emissivity is proportional to the 
number density of the stars, while the X-ray (bremsstrahlung) emissivity is proportional to the 
square of the hot gas n umbe r density, one obtains psir) oc PQ{r). This result was considered by 
[iVlathews & Brig hentil (^o3) as a general property of elliptical galaxies. Indeed, in the case of 
M 87 host galaxy, at the distances r > 18" one observes So(r) oc r^^-^ leading to ps(r) oc r 
and at the same time Sx(r) oc r^^-'^ leading to Pg(^) oc [Sx(r)/r]^/^ oc r 



-2.3 



-1.2 



(see above). This 



is in a good agreement with the expected behavior psir) oc Pcir). Therefore, hereafter we also 
assume that in the inner parts of M 87 the distribution of the hot thermal gas (i.e. of the pressure) 
follows closely the distribution of the stars (i.e. of the mass). We a lso assume for simplici ty a con- 



Di Matteo et al. 



stant temperature of the hot gas kTQ ~ 1 keV within r < 60" (see 
to the pressure profile simply proportional to the gas density profile, pcii") oc pcii"), i.e. 

-0.6 



2003h . leading 



Pair) = Po X 



1+ f: 



-0.6 



for r < Tb 
for r ^ Tb 



(1) 



with the normalization po = 1.5 x 10 ^ dyn cm ^ dYoung et al. l2002h . Here tb = 3" ^ 234 pc, and 



tk = 18" ~ 1.4 kpc. The resulting distribution of the gas pressure is shown in Figure 1. One can 
see that the gas pressure decreases from ~ 10^^ dyn cm~ ^ at r ~ 10 pc up to ~ 10~ ^° dyn cm~ ^ 



atr_~ 10 kpc. In addition, the pressure profiles adopted by 



Falle & Wilsonl(ll985h and 



Owen et al 



(119891) in their studies of the M 87 jet are shown for comparison. The former one is p{r) 
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HST-1 



D E F lABC 



E 

■D 



O 




-10 - 



2 3 

logr [pc] 



Figure 1. Profiles of the hot gas pressure in M 87 host galaxy, as evaluated bv lFalle&WilsorJ<1985l dashed line), IOwenetaljjl989l thin solid 
line), and in this paper (thick solid line). Circles indicate minimum pressure of the knots in the M 87 jet neglecting the relativistic correction 
(filled ones), and assuming the jet Doppler factor 5 = 2.7 (open ones). The circles disconnected from the others correspond to the HST-1 flaring 
region (the upstream edge of the HST-1 knot). In deprojecting distances between the knots and the active core we assumed the jet viewing angle of 
e = 20°. 

1.2 X 10"^ (1 + (r/re)^)"^/^ dyn cm-^ while the latter one is p(r) = 5.1 x 10"^° (r/kpc)-"-^^ 
dyn cm^^ between 0.7 and 2 kpc, and p{r) = 7.3 x 10^^° (r/kpc)~°'*^ dyn cm^^ for r > 2 kpc. 
Note that these two approximations imply a lower gas pressure than the one adopted by us. 

Figure 1 shows the de-projected positions of dif ferent knots (starting from HST-1 up to C) of 



Heinz & Begelman 



theM 87 jet^ assuming the jet viewing angle 9 = 20° dBicknell & BegelmaJl996[ 
199^, and the minimum pressure of these knots (a sum of the pressure due to ultrarelativistic ra 
diating electrons and due to the tangled magnetic field), 

/r V^W r 



Pmin Psq, e ~l~ ^B, eq 

(see 



3.4 X 10" 



Kataoka & Stawarz 



100 mJy 



0.3" 



-12/7 



5-^°/^ dyncm- 



(2) 



2005L and Appendix B). Here /r is the observed radio flux of the knot 



at 15 GHz, R is its observed knot's radius (assuming spherical geometry), and 6 is the knot's 

1 Knot HST-1: 0.8" - 1.2", knot D: 2.7" - 4", knot E: 5.7" - 6.2", knot F: 8.1" - 8.8", knot I: 10.5" - 11.5", knot A: 12.2" - 12.5", 
knot B: 14.1" - 14.5", knot C: 17.5" - 19". 
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Doppler factor. In Figure 1 we consider two jet's Doppler factors: 6 = 1 and S = 2.7. The later 
one is appropriate for the expected jet viewing angle ^ ~ 20° and the jet bulk Lorentz factor 



Heinz & Begelman 



1997h . We also took 



r ~ 3 — 5 on kpc scales (|Bicknell & Begel manlll996 : 
R = 0.3" (except the HST-1 flaring region, i.e., the upstream edge of HST-1 knot disconnected 
from the other knots in Figure 1, for which we assume R = 0.02") and used the knots' 15 GHz 



fluxes given by IPerlman et a 



flux of 3.8 mJy, as given in 



( 20011)^. For the fl aring region of HST-1 knot we took the 15 GHz 



Harris et al. 



(120031) for the quiescence epoch of this part of the jet. 
Note, that the first bright knot HST-1 is placed very close to re, i.e. the radius where the change in 
the ambient pressure profile between the central cusp and the unperturbed King-like distribution 
is expected to occur. In addition, downstream of this region, for tb < r < tk, the M 87 jet is 
overpressured in respect to the gaseous medium by a factor of a few, and even by more than an 
order of magnitude at the position of the brightest knot A further away. However, with the beaming 
effects included, the minimum pressure of the knots D, E, F and I is almost the same as the ambient 
medium pressure. Note also that the HST-1 flaring region is highly overpressured. 

An additional gaseous X-ray condensation in the center of M 87 host galaxy, linked to the 
observed in optical stellar cusp, increases a thermal pressure of the galactic medium (with respect 
to the 'pure' King-like profile) by as much as an order of magnitude at the distance ~ 10 pc from 
the core. On the other hand, a small volume occupied by this additional component implies only 
a small excess X-ray thermal luminosity. Figure 2 shows this light increase in the X-ray surface 
brightness profile. We calculate the X-ray surface brightness of host galaxy with and without this 
central component. Because the bremsstrahlung emissivity is proportional to the square of the 
thermal gas density, the appropriate surface brightness, as a function of the projected distance 
from the nucleus Vp, is 



Sx(rp)oc 1^ pl{^P + rl) dl , 



(3) 



where L 



rl and rx = 10 rc ~ 68.7 kpc (iLauer et al.ll992h . We assumed a constant 



temperature of the gaseous medium within the galaxy, and took the pressure profile as given by 
equation 1 with and without the central cusp. As shown in Figure 2, the additional central compo- 
nent increases only slightly the X-ray surface brightness, in particular by a factor of 2 — 3 within 
rp ^ 100 pc. 



Interestingly, 



Hardcastle et al 



— I 

( 2002) reported larger, by a factor of two, number of counts 



Knot HST-1 (total): 35.64 mJy, knot D: 161.54 mJy, knot E: 48.05 mJy, knot F: 144.9 mJy, knot I: 75.8 mJy, knot A: 1218 mJy, knot B: i 
mJy, knot C: 544.7 mJy. 
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from a thermal X-ray halo surrounding the central regions (rp < 500 pc) of FR I radio galaxy 
3C 31 than the number of counts expected from the pure (3 model fitted to the entire 3C 31 host 
galaxy profile. We believe that this excess can be related to a condensation of the hot interstellar 
medium in the central parts of the host galaxy, as discussed above. In the case of the M 87 source, 
however, it would be difficult to claim a presence of an analogous feature at r < 100 pc based 
on the available Chandra data, due to extremely complicated M 87 X-ray structure consisting 
of gaseous rings, voids, as well as due to strong non-thermal keV emission of the jet itself (see 



Feng et al.ll2004l:lDi Matteo et al.ll2003h . Because of such complexity, the thermal pressure profile 



proposed above should be really considered as a simple approximation only. For the purpose of 
the analysis presented below, it is however accurate enough. 

Let us note, that the spatial scale for the postulated here central gaseous condensation is very 
small when compared to the scale of the M 87 radio lobes. Thus, its presence does not contra- 
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diet with widely discussed disruption of cooling-flow atmospheres by kpc-scale radio outflows 



(see, e.g., 



Bicknell & Begelman 



199ft for the particular case of M 87 source). In fact, the sound- 



crossing time over the region with the spatial scale ~ tb is less than 1 Myr (for the interstellar 
medium parameters as co nsidered in this section), i.e . less than the lifetime of the inner lobes 
in the M 87 radio galaxy (IBicknell & Begelmanlll996l) . suggesting relatively short timescale for 
formation/regeneration of the central gaseous cusp. 



3 HST-1 KNOT AS A RECONFINEMENT SHOCK 



We do not intend to explain here the observed gradual coUimation of the M 87 jet in its inner- 
most parts. Instead we note that the initial coUimation of the b road nuclear outflow may be due to 
a dynamically dominating magnetic field dGracia et al.ll2005b . as the jets in active galactic nuclei 
are most likely launched by the magnetohydrodynamical processes. On the other hand, a dominant 
electromagnetic jet flux should be converted at some point to t he particle flux, sin ce the large-scale 



jets seem to be rather particle dominated (see a discussion in 



Sikora et a1. 



2005|) . Let us therefore 



speculate, that at sufficiently large distance from the nucleus — where the initial coUimation is 
completed — the relativistic jet in the M 87 radio galaxy is already particle dominated, and starts 
to expand freely. In such freely expanding jet, the pressure decreases very rapidly with the dis- 
tance, r, from the core. For example, in the case of cold jet matter the thermal pressure goes as 



Pj(rj oc r 



-27 



-10/3 foj. ^ ^ Jsanders 



1983r) . At the same time, the ambient gas pressure 



decreases much less rapidly: above we argue that in the M 87 galaxy one h as vair) oc r with 



Komissarov & Falle 



71 — .6 for r < 235 pc. Hence, as r] < 2, accordingly to the discussion in 

the initially free jet in M 87 certainly (i) will become reconfined at some point rg, (ii) 
will develop a reconfinement shock at its boundary, possibly leading to limb-brightenings of the 
reconfining outflow, and moreover (iii) the converging reconfinement shock will reach the jet axis 
at some further position along the jet, Vcr, beyond which the whole jet itself will come to a pressure 
equilibrium with the external gas medium. A simple evaluation of the reconfinement shock param- 
eters is presented in Appendix A for the cases of the jet m atter described by a non-rela tivistic equa- 



tion of stat e (hereafter 'co 



this context 



Sanders 



Komissarov 



1983 : 



d jet'), as done previously in 



Falle & WUsonlll985 



Wilson & FaUe 



198f 



Komissarov & Falld (|1997l) (see also in 

i 



Wilson 



19871: 



FaUe 



1991 



19941) . and also for an ultrarelativistic equation of state (hereafter ' hot jet') 



Reidetal 



( 19891) 



At what distance from the M 87 nucleus, r^, does the jet reconfinement start? 
noted that at the projected distance ~ 0.05" ^ A pc from the core the jet radio morphology 
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(opening angle, transverse intensity profile) changes. Further out, beyond ~ 0.1" ~ 8 pc, the 
jet brightness drops below the detection level, and then increases again at ~ 0.15" ~ 10 pc 
forming a weak radio knot L. Beyond this knot, the jet radio brightness decreases again, until 
~ 0.8" ~ 62 pc where a very bright knot, HST-1 appears (see section 1). Bearing in mind all 
the difficulties and uncertainties present in measurements regarding detailed morphology of the 
nuclear jet radio structure, we conclude that it is reasonable to identify rg with the jet region 
between 0.05" and 0.1" (i.e. 4 — 8 pc) from the M 87 center, and to assume that the jet thereby is 
already relativistic and particle dominated. Indeed, the initial — hydromagnetic by assumption — 
coUi mation of a broad nuclear o utflow seems to be already completed at smaller distances from the 



core. 



Giannios & Spruiti (120061) argue that the Poynting-flux dominated nuclear outflows in AGNs 



become kinetic flux dominated at distances > 10^ rg, i.e., in the case of M 87 radio galaxy, at 
about > 0.1 pc projected (for 6 = 20°), in agreement with our assumption. 

With ro ~ 0.05" — 0.1", one should expect the reconfinement shock to reach the jet axis at 
~ 3ro ~ 0.15" — 0.3" projected distance from the center in the case of a cold jet, or at 
Tcr ~ 10 ro ~ 0.5" — 1.0" in the case of a hot jet (see Appendix A). In other words, if the jet at 
ro is dynamically dominated by cold particles, r^r is expected to roughly coincide with the knot L, 
while for the ultrarelativistic jet matter — consisting of (mildly) relativistic particles plus magnetic 
field — rcr should rather be identified with the HST-1 complex. Again, noting all the rough approx- 
imations used by us to derive ro and rcr, below we argue that the latter interpretation is more likely. 
When the reconfinement shock reaches the jet axis, converging supersonic flow downstream of the 
reconfinement shock is expected to create the second stationary 'reflected' shock. This reflected 



shock is in turn diverging from the jet axis along the outflow, starting from rcr (IKomissarov & Falle 



19971) . The jet pressure immediately beyond rcr is therefore expected to be higher than the ambient 



medium pressure. This is qualitatively consistent with what is presented in Figure 1 for the HST- 
1 complex. Therefore we postulate that the extremely compact and overpressured HST-1 flaring 
point, present at the very beginning of the HST-1 complex, is placed at > rcr (thus favoring hot jet 
scenario), while the outer parts of the HST-1 complex — superluminal features characterized by 
the minimum pressure in rough equilibrium with the surrounding medium (see Figure 1) — can 
be identified with the region occupied by a diverging reflected shock further away from rcr. 

For a given rcr = '"cr.p/ sin6', where rcr,p = 0.8" = 62.4 pc is the projected distance of the 
HST-1 flaring region and 9 is the jet viewing angle, plus the adopted ambient pressure profile 
Pair) oc r^^ with r] = 0.6 and pq = 1.5 x 10^^ dyn cm^^, kinetic power of the jet implied by the 
model is 
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Figure 3. A total kinetic power of the M 87 jet as a function of the jet viewing angle d, implied by the model in which the reconfinement shock 
reaches the jet axis at the position of the HST-1 knot (solid line), and at the position of the knot L (dashed line). 



0.6 „1.4 



0.4 X 10^^ (sin 



ergs 



(4) 



(see Appendix A). The evaluated Lj is shown in Figure 3 for different jet viewing angles. Note that 
for 6 = 20° the im plied luminosity is L\ 10^'^ erg s~^, consiste nt with the jet power required 



to feed radio lobes (IBicknell & Begelman 



1996 : 



Owen et al 



20001) . Figure 3 presents also the jet 
kinetic power implied by the reconfinement shock position at r^^p = 0.15" = 11.7 pc (projected 
distance of the knot L), i.e. the location preferred in the cold jet scenario. In this case Lj is an order 
of magnitude lower. 

Previously, reconfinement shocks in FR I jets were proposed to be placed at much larger 
distances fror n the central engines, namely at the position of the brightest knots ~ 1 kpc from 



active nuclei dLaing & Bridle 



dFalle & WilsonI 



In fact, it is possible that beyond HST-1 complex (> 1.2") the M 87 jet 



2002h ^. In the case of the M 87 jet it would be then at knot A 



For a possibility of stationary reconfinement shocks in the small-scale jets of blazar sources see, e.g., 



Jorstad et al. 



2001). 
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breaks free again, and forms another reconfinement shock around 



10 X 1.2" = 12", i.e. ex- 



actly at the p osition of knot A. On t 



l e othe r hand, knots beyond HST-1 complex were successfully 



explained by iBicknell & BegelmanI (|l996|) as oblique shocks formed by helical modes of Kelvin- 



Helmholtz instabilities characterized by a g rowing amplitude alo ng the jet, disrupting finally the 



outflow near knot C (see also in this context 



Lobanov et al.ll2003h . For our following analysis, the 



discussion on the dynamics of the kpc-scale parts of the M 87 jet, beyond HST-1 is, however, not 
crucial. 



4 PHOTON FIELDS 



We evaluate energy densities of the ambient radiation fields along the jet axis, as measured in the 
rest frame of M 87 host galaxy at different distances from the center. First , we note that the op tical 



starlight emission is dominated by photons at frequencies of ~ 10 Hz (iMiiller et al 



20041), and 



that its emissivity profile is expected to follow the galactic mass (i.e., star) distribution. Hence, the 
emissivity is in a form 

\ -1.25 

for r < r-Q 



jstar(r) = jo X 



1+ 



-1.15 



for r ^ r-Q 



(5) 



with jo = 4 X 10^^^ erg s^^ cm^"^, corresponding to the /-band galaxy lumin osity (see section 2) . 
We integrate equation 5 along a ray and a solid angle with rx = 10^^ rc (see Stawarz et aklLoOS ) 
to obtain a profile of the starlight photons energy density for M 87, t/star = (1/c) / jstar [r) ds d^l, 
shown in Figure 4. Note that at distances r < 1 kpc it is roughly constant with < 10^^ erg cm^^. 
We can also evaluate the energy density of the X-ray photons for the observed X-ray emission of 
the hot gas with the temperature A; Tq ~ 1.5 keV in M 87. We note that the bremsstrahlung emis- 
sivity is simply proportional to the square of the gas number density, jism{r) oc riQ^r), and hence, 
with the assumed constant gas temperature, to the square of the gas pressure, jism('") oc PQ{r). 
By using the gas pressure profile given in the equation 1, and integrating jism(^) along a ray and 
a solid angle with the cluster termination radius ~ 1 Mpc, we obtain a distribution of the X-ray 
photons energy density, Uismii"), shown in Figure 4. Within the first kpc from the core the energy 
density of the thermal X-ray photons is higher than the energy density of the cosmic microwave 
background (CMB) photons, < 10^^^ erg cm^^. 

Energy density of the diffused radiation from stars and hot interstellar medium can be com- 
pared with the energy density of the synchrotron emission produced within each knot of the M 87 



Flaring HST-1 Knot inM87 Jet 1 5 




Figure 4. Energy density profiles for different radiation fields as measured in tlie rest frame of file M 87 liost galaxy by a stationary observer located 
at the jet axis. A thick solid line conesponds to the starlight emission, a dashed line to the thermal X-ray emission, a thin solid line to the CMB 
radiation, and a dotted line to the emission of the nuclear portion of the jet. Stars indicate energy densities of the internal synchrotron photons for 
different knots (neglecting relativistic corrections due to relativistic velocities of the emitting plasma). In deprojecting distances between the knots 
and the active core we assumed the jet viewing angle of 9 = 20° . 



jet. As discussed in section 1, synchrotron emission of the knots is peaked at optical frequencies. 
Thus, in order to evaluate the energy dens ity of the synchrotron photons within the jet, we take 



the optical fluxes measured at 10^^ Hz by 



Perlman et al 



( 20011) for aU the knots'^, except of the 



Harris etaf 



(|2003|) . corresponding to its 



HST-1 flaring point, for which we take 9 /iJy, as given in 
quiescence epoch. We also assume a spherical geometry for the emission regions, with the radius 
0.02" in the case of the knot HST-1 (considering only its flaring component) and 0.3" for the oth- 
ers. Figure 4 illustrates the resulting energy density of the synchrotron photons along the M 87 
jet (neglecting relativistic corrections), f/gyn = [^^o/o]/ {R^ c), where /o is the optical flux of a 
knot at z/q = 10^^ Hz. In deprojecting distances of the knots from the active core jet viewing angle 
6 = 20° was assumed for illustration. 



* Knot D: 59.5 /xJy, knot E: 16.2 ^jjy, knot F: 62.7 ^iJy, knot I: 28.6 ^iJy, knot A: 586 fj.]y, knot B: 306.8 ^ijy. knot C: 135.9 uly. 
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Finally, any stationary observer located at the jet axis is illuminated by the radiation produced 
within the active nucleus. Obviously, this emission is anisotropic, relativistically beamed into a 
narrow cone depending on the (unknown) bulk Lorentz factor of a nuclear jet. In this context 
we note, that the sub-pc-scale (r < 0.2 mas 0.015 pc) jet in M 87 has a slightly different 
posi tion angle than the la rge-scale jet in this source, with the misalignment in position angle ^ 
15° (IJunor & Birettalll995r) . For these reasons, it is not obvious that the large-scale jet is indeed 
illuminated from behind by the beamed nuclear emission and what is the beaming amplification 
of such an emission, i.e., if the jet flow at r > 0.015 pc from the center is placed within the 
beaming cone of the nuclear jet. If, however, jet misalignment can be neglected in this respect 
(because of relativistic and projection e ffects which make a pparent misalignment much larger 



than the real one) then, as discussed in 



Stawarz et al. 



( 20031) . energy density of the nuclear jet 

= Lnuc (2 Tnuc/^nuc) V(47r c), whcrc 



emission in the galactic rest frame along the jet axis is 
Lnuc is the synchrotron luminosity of the nuclear jet observed at some viewing angle 6, Fnuc is 
the bulk Lorentz factor of the nuclear jet and dnnc is the appropriate nuclear Doppler factor. For 
example, with 9 ~ 20° and Fnuc ~ 3 — 10 one obtains (2 Fnuc/^nuc)^ ~ 10 — 10^. In Figure 4, for 
illustration we assume the nuclear beaming correction factor ~ 10^, and note that the uncertainty 
in this approximation (fo r a fixed 6) can be more than ± one order of magnitude. We further take 



-^nuc = 3 X 10^^ erg s ^ (ITsvetanov et al. 



19981) characterizing steady state of the M 87 nucleus, 



obtaining thus a profile of f/nuc(^) shown in Figure 4. Let us mention, that Lnuc is peaked at the 
observed photon frequencies ~ 10^"^ — 10^^ Hz. 



In the rest frame of the jet, the energy densities of different radiation fields depend on the 
bulk Lorentz factor and inclination of some particular part of the jet. For example, the energy 
density of the starlight emission (as well as of the cluster and CMB photon field) are amplified 
in a plasma rest frame accordingly to oc F^. Relativistic corrections also decrease the comoving 
energy density of the synchrotron radiation accordingly to oc (as appropriate for a stationary 
shock feature). Finally, the nuclea r emission in the rest frame of the outer jet is decreased by a 
factor (2 r)~^ ( Stawarz et akllOOB ). Note, that even varying the jet viewing angle alone influences 
deprojected distances of the jet features and therefore the energy densities of the galactic and nu- 
clear radiation fields. In a framework of our model, HST-1 flaring point corresponds to a compact 
region just downstream of the reconfinement/reflect shocks system. Estimation of the appropriate 
bulk Lorentz factor of the radiating plasma is not trivial in this case, because we need to consider 
the oblique shock geometry. In particular, the jet matter downstream of the reconfinement and 
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Figure 5. Energy densities of different photon fields measured in a rest frame of tire HST- 1 flaring region, as functions of the bulk Lorentz factor of 
this part of the jet, F, for different jet viewing angles. Thick solid lines correspond to the starlight emission. Dashed lines correspond to the internal 
synchrotron emission of the knot, for the jet inclination 9 = 10°, 20°, 30°, and 40° (from bottom to top, respectively). Shaded region indicates 
energy density of the equipartition magnetic field for 6 = 10° — 20° . Dotted lines illustrate comoving energy density of the nuclear photons for 
e = 40°, Tnuc = 10 (upper curve), and 6» = 10°, T^uc = 3 (lower curve). 

reflected shock fronts may still be relativistic, depending on the distance from the jet axis (see 
Appendix A). 

If the HST-1 flaring point corresponds to a significantly decelerated portion of the jet matter 
located at the very jet axis close to Vcr, the comoving energy densities of the starlight and syn- 
chrotron photon fields are comparable, as presented in Figures 4 and 5 for different jet viewing 
angles. With an increasing bulk Lorenz factor the energy density of the starlight emission increases 
from about ~ 3 x 10~^° erg cm~^ for F ~ 1 up to ~ 3 x 10~^ erg cm~^ for F ~ 10 (Figure 5). A 
shift in the deprojected position of the HST-1 flaring region due to a different jet viewing angle is 
of negligible importance, because the entire HST-1 complex is located within the central plateau 
of Ustarir). The internal energy density of the jet synchrotron emission initially decreases with a 
growing velocity of the emitting region, but for the larger values of F it increases again, since for a 
given 9 ^ 10° the appropriate jet Doppler factor 5 decreases with an increasing (large) F. One can 
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however conclude, that for a wide range of parameters shown in Figure 5 (namely, 6 = 10° — 40° 
and r 1 — 10), in the rest frame of the HST-1 flaring region, the starlight emission is expected to 
dominate over the internal synchrotron photon field. On the other hand, one should be aware that 
as the HST-1 flaring point is unresolved (while at the same time synchrotron energy density goes 
with the emission region size as oc R~'^), the estimated U'^^^ should be considered as a lower limit 
only. In addition, here we only considered a quiescent epoch of M 87. The biggest uncertainties 
correspond however to a photon field of a nuclear jet illuminating HST-1 knot from behind. Figure 
5 illustrates two extreme cases for f/^^^ at the position of this knot, corresponding to the nuclear 
Lorentz factor Fnuc = 10 and 9 — 40°, and also to Fnuc — 3 and 9 — 10°. The estimated energy 
density of the synchrotron emission of the nuclear jet varies by a few orders of magnitude (!) for 
these two examples, and may exceed or be much smaller than the other components, C/^g^ and 
[/gyjj. Figure 5 shows also for a comparison the energy density of the equipartition magnetic field, 
[/b = {Bl JM) 5"^°/^ with = 10-3 G, for ^ = 10° - 40°. 



5 HST-1 KNOT AS A TEV SOURCE? 

Let us suppose that the active core of M 87 experienced at some moment an outburst, resulting in 
the flare of its synchrotron emission and ejection of a portion of the "jet matter" with the excess 
kinetic power (when compared to the steady-state epoch of the jet activity). Both photons and par- 
ticles travel along the jet, arriving at some time to the location of the HST-1 knot, where the recon- 
finement shock formed within a steady jet reaches the jet axis. Flare synchrotron photons emitted 
by the active nucleus are then comptonized to TeV energies (hereafter TC/nuc' process), while 
the excess jet matter shocked around causes synchrotron (hereafter 'SYN') and the additional 
inverse-Compton brightening of the HST-1 flaring region. As discussed in the previous section and 
also below, this additional inverse-Compton brightening should be dominated by Compton scat- 
tering of the starlight emission (hereafter TC/star' process) or synchrotron self -Compton process 
('SSC'). We note, that some short time-scale variations of the emission coming from the nucleus 
in its high state can be imprinted in the observed IC/nuc flux of the outer parts of the jet. In ad- 
dition, due to different velocities of the nuclear photons and particles, TeV flare resulting from 
the IC/nuc process in HST-1 knot is expected to lead SYN, IC/star and SSC brightening of the 
HST-1 flaring region by some time At. Note that the increase in the seed photon energy density 
for a given synchrotron flux (i.e. for a given particle energy density) results in an increase of the 
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inverse-Compton flux only if the electrons involved in the inverse-Compton scattering are weakly 
cooled by radiative losses ('slow cooling regime'). 

Assuming that the observed sub-luminal velocities of the jet features observed between the 
core and the HST-1 knot are only pattern velo cities, not reflecting the true bulk velocity of the jet 
spine (see a discussion in Dodson et alJLoOS ). and that this true bulk velocity is highly relativistic 
/3nuc = (1 — ^nuc)^^^^ ~ 1, the appropriate delay time difference is roughly 

At--^-^-— ;5L^100(sin^^r^r- yr, (6) 

where = r sinO = 62 A pc is a projected distance of the HST-1 flaring region from the core. 
For example, period At ~ 6 yr between presumable maximum of the TeV emission (1998/1999) 
and the observed maximum of the synchrotron emission of the HST-1 knot (2005) is consistent 
with the jet viewing angle 9 ~ 10° for Fnuc ~ 10, with 9 ~ 20° for Fnuc ~ 7, and finally with 
9 ~ 30° for Fnuc ~ 6. The assumed hypothetical nuclear flare should be observed some tfl years 
before the IC/nuc flare of the HST-1 knot, where 

T 1 — POS 9 

tfi ~ - - — ^ ~ 200 (1 - cos 9) (sin 9)-' yr . (7) 
c smfc* 

For example, tfl ~ 20 yr for 9 ~ 10°, t f, ~ 35 yr for 9 ~ 20°, and ~ 55 yr for 9 ~ 30°. 
We note in this context, that interestingly De YoungI (1197 11) reported radio flare of M 87 nucleus 
in 1969-1971. If — again for illustration — one identifies the considered nuclear flare with this 
event, then equations 6 and 7 imply ~ 16° and Fnuc ~ 8. 

In the rest frame of the HST-1 knot, assuming moderate bulk velocity and jet viewing angle, 
the energy densities of the starlight, nuclear and internal synchrotron photons are peaked at similar 
photon frequencies 10^'' — 10^^ Hz. Thus, electrons upscattering all these photons to the observed 
TeV energies are mostly the slowly cooled ones, with energies ~ 10^ mgC^ (the optical spectral 



Perlman et al 



2003h . The resulting TeV fluxes 



index of HST-1 knot is consistent with «o ~ 0.6; 
due to the IC/star, IC/nuc and SSC processes are then produced in the transition between Thomson 
and Klein-Nishina regimes. Hence, for a rough evaluations one can approximate the expected ob- 
served TeV fluxes by Ljc/seed ~ /±,iso (^seed/f^i?) ^o, where Lq = 47rrfL [t'o/o] is the observed 
optical luminosity of the HST-1 knot, while ^Tg'eed the comoving energy densities of the 

appropriate seed photons and the equipartition magnetic field. Function = /±, i^n(T,9) arises 



due to possible anisotropy of the seed photons in the emitting region rest frame (see 



Stawarz et al 



20031 and Appendix B). Using the observed radio fluxes /r of the HST-1 flaring region as mea- 
sured at 15 GHz, the optical flux /o at 10^^ Hz, and the emitting region size R in arcseconds, one 
obtains (see Appendix B) 
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Figure 6. Expected TeV emission of the HST-1 flaring region in 1998, due to IC/nuc (dotted line), IC/star (solid line), and SSC (dashed lie) 
processes, as a function of the bulk Lorentz factor T of this part of the jet assuming 9 = 20° . Shaded regions indicate the appropriate luminosity 
ranges for 6» = 20° - 30°. 

. 7.3 X 10- (A) (_|) (_^) e.g s- . (8) 



and 



. ,4 X 10- (A) (_^)''" 

X (sin^)2(l-cos^)^ ergs"^ (10) 

In the above we assumed that at every moment (i.e., for a given /r and /o) the HST- 1 flaring region 
is in equipartition regarding energies of the ultrarelativistic electrons and the magnetic field. We 
also introduced internal nuclear luminosity Lg, which should correspond to the assumed nuclear 
outburst, and not to the steady-state discussed in the previous section. 

Different constraints presented and cited in previous sections suggest the most likely jet view- 
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Figure 7. Expected TeV emission of the HST-1 flaring region in 2004, due to IC/nuc (dotted line), IC/star (solid line) and SSC (dashed line) 
processes, as function of the bulk Lorentz factor F of this part of the jet assuming d = 20°. Shaded regions indicate the appropriate luminosity 
ranges for 6» = 20° - 30°. 



ing angle of ~ 20° . In Figure 6 we present the expected TeV emission of the HST- 1 flaring region 
in 1998 assuming = 20°, and resulting from IC/nuc, IC/star and SSC processes (equations 8- 
10), as functions of the bulk Lorentz factor of this part of the jet. For illustration, shaded regions 
indicate also the appropriate lu minosity expected for 9 = 20° — 30°. Here we took /r = 3.8 mJy, 



20031). The expected TeV IC/star emission is in this 



/o = 9 /iJy and R = 0.02" (iHarris et al. 
case Lic/star < 10'^^ erg s~^, and the corresponding SSC emission is even lower. However, the 
IC/nuc emission could eventually account for the HEGRA detection (L^ > 10^° erg s^^) only if 
the bulk Lorentz factor of the HST-1 flaring region was F ~ 2 and the assumed nuclear flare was 
characterized by F^^^Lg > 3 x 10^® erg s~^ (as taken in Figure 6). We note, that for such param- 
eters the comoving energy density of the nuclear photons, U[^^^ > 10" erg cm~^, dominates over 
the comoving energy densities of the magnetic field and the starlight emission {U!q ~ 10"^ erg 
cm"'^ and U'^^^^ ~ 10~^ erg cm"'^, respectively). Hence, cooling of the considered ~ TeV energy 
electrons is mainly due to the IC/nuc process. We also note, that the obtained above value of the 
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preferred bulk Lorentz factor F ~ 2 — 3 refers in our model strictly to the compact (unresolved) 
and decelerated portion of the outflow placed at the jet axis immediately after Td-, which is respon- 
sible for production of the flaring emission, and not to whole outflow at the same distance from the 
nucleus. In fact, oblique geometry of the reconfinement shock implies that the other parts of the 
jet, located further from the jet axis, may suffer much less deceleration, and thus that the average 
bulk Lorentz factor of the whole outflow may be higher than the one characterizing HST-1 flaring 
region. 

With the preferred ^ ~ 20° - 30° and T ~ 2 - 3, the Doppler factor of the HST-1 flaring 

region is 5 ~ 2 — 3, while the jet-counterjet synchrotron brightness asymmetry is /j//cj ~ [(1 + 

(3 cos 6)/{\ — [3 cos Oyi^-^ ~ 10^ — 10^ (still constistently with the observational limits). Assuming 

Tnuc ~ 10 (in agre ement with values usually derived by means of modeling broad-band emission 

of BL Lac objects; Urrv and PadovanJl995 ). the Doppler factor of the nuclear M 87 jet would be 

roughly 5nnc ~ 1- This implies that the nuclear outburst assumed in our model would be observed 

with the isotropic luminosity ~ (5nuc/rnuc)^ (rnuc-^fi) ~ 3 x 10^^ e rg s~^. This is higher t han 

the observed luminosity of the M 87 nucleus in its steady-state epoch (|Tsvetanov et al.lll998|) by 

only a factor ~ 10. We note, that order-of-magnitude flares on time scales of years are common in 

blazar sources. In addition, as mentioned above, synchrotron e mission of the HST -1 flaring region 

J '. II 

has increased between 1998 and 2005 by a similar factor ~ 50 (IHarris et al.ll20 06). This is another 
indication that the model presented here is self-consistent (if only roughly L^mc oc Lj). Moreover, 
we expect characteristic timescale for the variability of the emission produced within the HST-1 
flaring region tvar ~ R/c5 < 1 yr for i? < 0.02", again in a rough agreement with the observed 
one at radio, optical. X-ray, and 7-ray frequencies. 

Let us investigate next the expected TeV emission at some later time, during the synchrotron 
flare of HST-1, when the nuclear seed photon energy density in the knot's rest frame has decreased 
significantly. Figure 7 shows the expected TeV IC/nuc, IC/star and SSC luminosities again for 
^ = 20° - 30° and R = 0.02", but this time with /r = 4 mJy, fn = 200 u Jy taken to iUustrate 
synchrotron continuum of HST-1 flaring region in 2004 (IHarris et al.ll2006l) . and Tl^^L'f^ = 10^^ 
erg s~^ characterizing the quiescence nuclear emission (factor of 30 below the high activity epoch 
considered above). In such a case, for F ~ 2 — 3 one expects Lic/star ~ (1^6) x 10^^ erg s^^ and 
-^ssc < 10^^ erg s"^. Also Lic/nuc < -Z^ic/star except for the small bulk Lorentz factors (F < 2) 
and large jet viewing angle (6 ~ 30°). The considered parameters imply now ~ 10^^ erg cm^^ 
and U^^^ ~ 3 X 10~^ erg cm^^, i.e. that radiative cooling of the TeV energy electrons is mainly due 
to their synchrotron emission. We finally note that the synchrotron emission of the HST-1 flaring 
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region has increased between 2004 and 2005 by a factor of 2 — 3. Thus, in a framework of the 
presented model, we expect also the TeV flux due to the IC/star process to increase in 2005 when 
compared with the 2004 level. 

6 SUMMARY AND CONCLUSIONS 

Supermassive black holes present in centers of active galaxies are known to influence trajecto- 
ries of nearby stars, and to create in this way central stellar cusps observed by Hubble. Here we 
propose that the distribution of the hot gas within ellipticals follows closely distribution of the 
stars not only in the outer parts of the galaxies, as observed in a number of such systems, but 
also in the innermost parts. If this is the case, then one should expect excess of thermal pressure 
(when compared to the pure /3-type profile of the gas number density) within ~ 100 pc from the 
galactic center. The resulting small excess in X-ray surface brightness due to free-free radiation of 
the hot gas seems to be required to explain some Chandra observations. This additional gaseous 
component can also result in a stronger confinement of the jets, leading to formation of stationary 
reconfinement/reflected shocks within the outflows. We propose that in the case of the M 87 radio 
galaxy, HST-1 knot present at ~ 100 pc from the center can be identified with the downstream 
region of such a reconfinement/reflected shock. In particular, we argue that stationary, compact 
(i? ^ 1 pc), variable (on the time scale of, at least, months and years), and overpressured (by a 
factor > 10) HST-1 flaring region is placed immediately downstream of the point where the con- 
verging reconfinement shock reaches the jet axis ('reconfinement nozzle'). Thereby some portion 
of the hot relativistic jet decelerates from highly relativistic to mildly relativistic bulk velocities 
(from bulk Lorentz factor F ~ 10 down to F ~ 2 — 3), while other parts of the jet (placed further 
away from the jet axis) are expected to decelerate less strongly due to a larger angle between the 
upstream bulk velocity vector and the shock normal. The liberated bulk kinetic energy of the out- 
flow is transformed at the shock front to the turbulent magnetic field energy (consistently with the 
decrease in the degree of linear polarization observed in HST-1 knot), and, in similar amount (by 
assumption), to the ultrarelativistic particles. 

Although the reconfinement/reflected shock structure is stationary in the observer's rest frame, 
variations and changes in the central engine lead inevitably to flaring of this part of the outflow, in 
particular when the excess particles and photons emitted by the active nucleus in its high-activity 
epoch and traveling down the jet arrive after some time to the reconfinement nozzle. In a frame- 
work of this scenario, one should expect firstly high-energy 7-ray flare due to comptonization of 
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the photons from the nuclear outburst, and then, after some delay depending on the bulk velocity 
of the nuclear jet, synchrotron flare due to excess nuclear particles shocked at the nozzle. This de- 
layed synchrotron flare could be accompanied by the subsequent inverse-Compton brightening due 
to upscattering of the ambient radiation fields by the increased population of the ultrarelativistic 
particles. It is tempting to speculate that such a sequence of events was in fact observed in HST-1 
flaring region, especially as for a realistic set of the jet parameters the evaluated radiative fluxes 
are in agreement with the multiwavelength observations performed between 1998 and 2005. 
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APPENDIX A: RECONFINEMENT SHOCK 



In the rest frame of a shock, the relativistic shock jump conditions can be written as 



(Al) 



w_ (3. = w+ rl /?+ , 



and 



n_ r_ = ri+ r+ 



(A2) 



(A3) 



where velocities /3_ and (3^ refer to the normal co mponents of the upstre am ('— ') and downstream 

.,|Kii 



('+') bulk velocity vectors, respectively (see, e.g.. lKirk & Duffylll999|) . Here w is the proper en- 
thalpy of the fluid, p is its isotropic pressure, and n its the proper number density. Let us consider 
first the case of the upstream cold plasma dominated dynamically by the rest energy of the particles 
with a mass m, with the thermal pressure negligible, i.e. the enthalpy being approximately equal 
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to the proper rest-mass energy density /i_ = mn_ c^, namely W- = fi- + 7P-/ (7 — 1) ~ 
where 7 is the appropriate adiabatic index. One can find that in such a case 

P+ = /._r2_/3! (l-^). (A4) 



Now let us consider a supresonic jet which breaks free at some distance from the central engine, 
an d next experiences reconfine ment by the ambient medium starting from the distance tq. Follow- 



ing 



Komissarov & Falld (Il997h . we denote by the angle between the tangent to the converging 



reconfinement shock at some given distance r > ro, and by the angle between the pre-shock 
jet bulk velocity vector close to the shock at the same distance r. Note, that by the definition 
tanip = —dz/dr and tan0 = z{r)/r, where z is the distance of the reconfinement shock from 
the jet axis at given r. We also assume that both angles are small, i.e. tan^ ^ ip and tan0 ^ (p. 
As the reconfinement shock is stationary in the observer rest frame, one has 

/3_ = /?j sin + 0) , (A5) 

where (3^ is the pre-shock jet bulk velocity, and Fj = (1 — fi^"^)^^/"^ = r_ is the pre-shock 
bulk Lorentz factor, and, obviously, /xj = The jet luminosity is Lj = WjT^ Pjcn R"^ ^ 
/ij rj/?j cvrr^ tan^ $, where = r tan$ is the radius of the free jet and $ is the pre-shock 
(free) jet opening angle. Taking the external pressure of the ambient gaseous matter pair) = 
Po {r/r-Q)~'^, by means of the condition p+{r) = pcir), one obtains an equation for the distance 
of the reconfinement shock from the jet axis 

^ = f_tt^(2-.)/2^ (A6) 

dr r 
where 

^^fPolhl^^Y' (A7) 



and we expressed the term (1 — as a parameter Ci. With the initial condition z{ro) = zq = 

ro tan $, the solution to the above equation, 

z{r) = r tan $ - ^ r (^r^^~v)/2 _ ^{2-'y)/2 j ^ (A8) 
implies that the reconfinement shock reaches the jet axis at the distance 



4 Po'^bCtt 
Note, that at tq one has 



(A9) 
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where M.-^ is the relativistic Mach number of a free jet, and /?s,j = (1 — Tg*^)^/^ is the jet sound 
speed in c units. For a cold jet matter considered here Fg j /5sj ~ /5sj = (7Pj//ij)^/^. Since at tq 
jet pressure equals external gas pressure, ]5j(ro) = Pg{^o)^ one can therefore find that 



47 ,7 2-r, 



^ (7-1: 

This, together with the equation for r^, gives the condition 



n l/{2-r?) 



(7 - 1)^ 

With ^ = 5/3 one obtains r^,/rQ = [2.625 (2 - r])'^]^^^^'"'^ 



(All) 



(A12) 



for ~ 0.7 (see 



Komissarov & Falle 



1997h . 

Now let us consider an analogous case as before, but with a jet matter described by an ultrarel- 
ativistic equation of state, w = Ap (i.e., with 7 = 4/3). In this case the upstream pressure cannot 



be neglected anymore, and by combining equations A1-A2 one obtains for a relativistic jet 



(A13) 



where C2 = ^ ~ iP+/ P-) = 1 — (1/3/9^) ^ 0.65. With the appropriate expression for the jet 
kinetic luminosity, Lj = 4pj F| /5j c vr tan^ $, one obtains again 



J 

l/(2-»?) 



, (A14) 

In addition, in the case of the ultrarelativistic equation of state the sound speed is /3c, j = l/"\/3, 
and hence by means of expression A 10 in a form 

3^2 



tan$ 



(A15) 



one obtains condition rcr/ro = [18 (2 — r]Y C2] 



l/(2-r?) 



2ll/(2-ry) 



APPENDIX B: RADIATIVE FORMULAE 



For a given radio flux /r as measured at some observed radio frequency z/r, and for the observed 
emission region size R, the intensity of the equipartition magnetic field evaluated ignoring rela- 
tivistic correction is 



-Beq,5=l OC 



(y. £ ( — f 



Q+1/2 



LX / 



R- 



-|2/7 



(Bl) 



where u^^^ and z/max are the minimum and maximum frequencies of the synchro tron continuum . 



assumed to be a simple power-law characterized by a spectral index a (see, e.g.. 



Longaii 



1994) 
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As discussed by 
taking z/max > t'l 



Stawarz et al. 



( 2003r) . relativistic corrections give B, 



cq 



~ z/R and a > 0.5, one obtains 



Ben 8=1 5 Hence, 



(B2) 



Kataoka & Stawarz 



(see in this context 
energy density t/^ o c B\ oc f^^ i ? ;^^/^ ^ 



20051). This gives the comoving minimum magnetic field 



As discussed in 



Stawarz et al 



(|2003r) . the observed inverse-Compton luminosities (produced 



in the Thomson regime) can be simply evaluated as 



^IC/seed 



it, iso 



TJI 

seed T 



(B3) 



where Lgyn is the observed synchrotron luminosity, f/ggcd the comoving energy densities of the 
seed photons, while /±,iso = /±,iso(r, d) is the function of the kinematic jet parameters arising 
due to possible anisotropy of the seed photon fields in the jet rest frame. In section 5, the observed 
bolometric synchrotron luminosity is approximated by the optical one, Lsyn oc z^o/o- In section 5 
we also assumed that at every moment (i.e., for a given synchrotron flux), the emission region is 
in the equipartition regarding energies of the radiating electrons and the magnetic field. 



In the case of the synchrotron self-Compton emission, /iso = 1 and U!.^^ oc /o -R '^6 ^, leading 



to 



T ^ 4/7 r2 p-2/7 r-11/7 



(B4) 



If the comptonisation of the starlight emission is considered, /+ ~ (S/T)'^ and U^^^^ oc F^. Hence, 



(B5) 



Finally, for the comptonisation of the nuclear emission illuminating the jet from behind, we have 
/_ ~ 52^2 (1 _ pQg^)2 and U;^^^ oc (L^F2^jF-2 (sin^)2, where the factor (sin^)^ is due to 



deprojecting the observed distance of the emission region. This leads to 

f-'/' fo R'^'' (4FL) 5'"' (sin ef (1 - cos Of . 



(B6) 



LlC/mic OC 

The approximations derived above allow us to estimate in a simple way the expected 7-ray 
fluxes for a given /r, /q, R, and the kinematic parameters of the jet. 
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